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Study guide 


Before reading this unit you should read those sections of Unit 25 which refer to 
the calculation of partial derivatives. 


You should read Sections 1-3 before watching the television programme, but you 
can choose to read Section 4 either before or after the programme. 


1 Introduction and definitions 


1.1 Introduction 


This unit is devoted almost entirely to what can truly be described as one of the 
‘gems’ of late eighteenth and early nineteenth century mathematics—the solution 
of partial differential equations with given boundary conditions—although in the 
space available we shall be able to examine only one simple case in any detail. 


My objective will be to give you a glimpse of the power of the methods which you 
have studied so far—and incidentally, and more prosaically, to enable you to 
revise some of the topics of this course in what I hope is an interesting fashion. 


The subject which we shall study arose in a period when a number of great 
mathematicians were turning their minds to the investigation of various questions 
arising in what we should call applied mathematics. To their credit, the men of this 
age would have drawn no clear distinction between what was ‘pure’ and what 
‘applied’; they were simply intent on obtaining answers to interesting questions, 
and in searching for these answers they incidentally laid the foundation for 
modern mathematical analysis as well as for modern physics. 


A catalogue of the names associated with the subject reads like an honours list: 


Euler, Laplace, Lagrange, Poisson, Cauchy, the Bernoullis, Fourier, and many See, for example, Chapter 22 
of Mathematical Thought from 


others. Before we become involved in the details of the single example which we 

ee Ancient to Modern Times, by 
have time to study, it will perhaps give you a deeper understanding of the Morris Kline (Oxford 
relevance of the topic if I say just a little about the range of subjects which these University Press, 1972). 


men considered. 


A number of apparently distinct areas of. physics — electrostatics, magnetism, 
hydrodynamics, aerodynamics, heat conduction and wave motion— use models 
which involve differential equations with partial derivatives. Such an equation is 
called a partial differential equation. The equation 


a0 _ 00 
6x? ~—sst 


is an example of a partial differential equation, and it can be used to model the 
flow of heat along a thin metal bar; in this case @(x, t) is the temperature at 
position x at time t. Having deduced that 6 satisfies such an equation, we are very 
little further forward unless we are able to determine 6 explicitly, using any 
available information about the temperature at some particular point in time and 
how heat may be lost from the system. 


The break-through in such problems came with the development of trigonometric 

series and the publication of Joseph Fourier’s Theorie Analytique de la Chaleur 

(1822), which is the source of many important methods in mathematical physics, in 

particular for problems involving partial differential equations under given 

boundary conditions. Certainly the subject had been advanced by Euler, 

D’Alembert, David Bernoulli and others, but it was Fourier who made the 

situation absolutely clear by establishing the fact that an ‘arbitrary’ function (in 

fact a piecewise continuous function) could be expressed in the form Unit 31, Section 2 


>> (A, cos nax + B, sin nax). 
— n=0 
Despite the support of Euler and Bernoulli, the idea was so original when first 
presented in 1807 that it raised some vigorous opposition, notably, it is said, from 
Lagrange (the leading French mathematician of the day). The reluctance of some 
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mathematicians to accept Fourier’s results must certainly have been due in part to 
the fact that it called into question the very meaning of the fundamental concept 
of a function. Ultimately, the ‘Fourier series’ was universally adopted, and the 
questions which it raised spurred mathematicians to re-examine the foundations of 
the calculus and to restore it in its modern form as mathematical analysis 


The essential elements of Fourier’s method are easy to describe. First derive the 
appropriate partial differential equation, then assume that the required function 
can be expressed in terms of a Fourier series, and finally use the equation and the 
initial conditions to determine the constants in the series. This is precisely the 
outline of the core of this unit; but, before we start on this programme, I should 
like to look at the problem that we shall attempt to solve in more detail. 


Imagine a coiled spring stretched tightly between two fixed points A and B 
(Figure 1). Now suppose that a point in the middle of the spring is pulled to one 
side, towards B, and then released, while the ends are kept fixed (Figure 2). The 
spring will start to vibrate, with each point of the spring moving to and fro in 
some fashion. Our task is to determine how each point will move, and to discover 
where it will be at each subsequent instant of time. In our solution of this problem 
there will be three main factors: 


(i) the physical properties of the spring, which will give rise to an equation 
known as the wave equation; 


(ii) the physical constraints applied to the spring, namely that it is fixed at 
both ends, which give rise to the so-called boundary conditions; 


(111) the initial behaviour of the spring, in other words, what we do to the 
spring in order to make it move. This initial behaviour produces the initial 
conditions for the problem. 


The exact meaning of these terms will become clear as we proceed through the 
unit. 7 


Although I have chosen to focus on the motion of a vibrating spring, the methods 
which we shall investigate are very general, and apply to the equations used for a 
wide variety of physical situations including heat flow, the sound produced by 
musical instruments and the motion of fluids. 


The essence of each of these physical situations is that it can be modelled by a 
partial differential equation, and so I must explain in more detail what I mean by 
such an equation. 


1.2 Partial differential equations 


Previously in this course we have solved a differential equation such as 


ey dy 

— + b— = 4), 

<3 ve = ay + cy 

which is an example of an ordinary differential equation, and is characterized as 
such by the fact that the unknown function y(x) is a function of the single variable 
x. We are now intending to discuss differential equations which apply to functions 
of two variables. Such a function, U(x, t) say, can vary both because of changes in 


x and because of changes in t. 


: ee ee) een : : 
As an illustration, consider the equation x —- = t—. The function U(x, t) = e* is 


: Ox Ot 
a solution of this equation, as we can easily verify. We have 
as = te" - and it =o te, 
Ox ot 


and it follows that md aie phe and so 
Ox Ot 


U(x) = & 


A AWAWWYVWVW/VIIIVIW--_ B 


Figure 1. A coiled, stretched 
spring. 


A SASAILAIA9Wv-- B 
aD 

Figure 2. The same spring 

with the middle pulled 

towards B. 


We shall see later how this use 
of the term ‘boundary 
conditions’ relates to the use 
in Unit 6, Section 3. 
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is a solution of the partial differential equation 


0U oU 
sae ot 
As you may have guessed, the general problem which we should like to solve is to 
discover all of the functions which satisfy a given partial differential equation, and 
this is the process which we call solving the equation. Generally speaking, solving a 
partial differential equation is by no means easy, and the following illustration 
may help you to understand why. 


Suppose that we are given the rather simple looking partial differential equation 
CU, i). 
eee 


The left-hand side is the partial derivative of U(x, t) with respect to x keeping t 
constant, so to find U we may keep t constant and integrate with respect to x. We 


then obtain 
U(x, t)=x+C. Cie 


However, the constant C is only constant if t is kept fixed, and for different values 
of t the value of C will be different, so that C is a function of t. Thus the most 
general solution is 


U(x,t)=x+ f(t) 


where f is an arbitrary function. You can verify that any such function U(x, f) is a 
solution by differentiating partially with respect to x. 


:. 


t no longer 
fixed 


So you seé that while the solution of an ordinary differential equation will involve 
arbitrary constants, the solution of a partial differential equation will involve 
arbitrary functions. 


As we increase the order of the partial derivatives in a partial differential equation, 
so we introduce more arbitrary functions into the solution. Thus the partial 
differential equation 

6*U 

Oxot 
might be expected to contain two arbitrary functions in its solution, and the 
following argument shows. that this is indeed so. 


A 


oU 
ae 
arbitrary function (of one variable). Now integrating this equation with respect to 
t (for constant x), we get 


Integrating with respect to x (for constant t) we see that f(t), where f is an 


U(x,t) = | f(tdt + g(x), 


where g is another arbitrary function. This solution can be written in the form 
U(x, t) = h(t) + g(x) 
where both h and g are arbitrary functions. 


In order to make any progress it will be necessary for me to restrict the class of 
partial differential equations to be discussed. Partial differential equations are 
classified in much the same way as ordinary differential equations (see Unit 6), so 
that a linear second-order partial differential equation in two independent variables 
will be one of the form 


6?U 67U 6*U oU oU 
a(x, er + b(x, ar + c(x, ioe + d(x, a + e(x, 1) + f(x, t)U = g(x, t) 


where a, b, c, d, e, f and g are known functions of x and t, and at least one of a, b 
and c is non-zero. 
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This is still far too wide a class to tackle in one unit, and in fact our concern will 
be with equations which (in addition to obeying the above conditions) are 
homogeneous (i.e. g(x, t) is the zero function) and constant-coefficient (i.e. a, b, c, d, 
e and f are all constants). Thus, our attention will be confined to equations of the 
form 

6?U 6°U 6*U oU oU 


Oa aa oe te Fi (1) 


where a, b, c, d, e and f are constants. 


These are the analogues of the ordinary differential equations of the form 


d°y. dy 
ee ae 
a2 + x +cy=0 


which are discussed in Section 1 of Unit 6. Theorem 1 of that section also has the 
following analogue for partial differential equations, which is of extreme 
importance to us for our work later in this unit. 


Theorem 1 


If u and v are two solutions of a linear homogeneous partial differential equation, 
then for all values of A and B, the function Au + Bv is also a solution of the 
equation. 


Proof (for the second-order constant-coefficient case) 


Let U = Au + Bo, where u and v are solutions of Equation (1). Then 


07U 0? 
es —— 5 Au + Bv) 
zy 07 v 
= aA a + ab 3 
07u 07v 
A as] + 2 as] 


-and similarly for the other terms in the equation. Thus 


a hu et 
eee fy 
eo t's oe eG tl 


67u 6*v 
= Aas ++ ful + Bla + : + fi 

=O (since u and v are solutions of Equation (1)), 
so that U is a solution of Equation (1). | 


In fact this theorem holds (as stated) for all linear homogeneous partial differential 
equations; the general proof is longer than the above but very similar in its 3 
method. | 


Before we go any further, I must make sure that you are confident with partial 
derivatives, so if you have any doubts about your ability try the following 
exercises. | 


Exercise 1 


OF OF 
Find = and oy for each of the following functions F(x, y). 


(i) F(x,y)=— 

y 
(ii) F(x,y) = Ax? + 2Bxy + Cy? 
(iii) F(x, y) = sin(2x + y) 


(iv) F(x, y) = @(2x + y), where @ is an arbitrary differentiable function of one variable. (Use 
the notation @¢’ for the derived function of @.) 


(v) F(x, y) = b(ax? + by?) where ¢ is an arbitrary differentiable function of one variable. 


MST204 32.1 


Exercise 2 3 
Write down simple expressions for the following: 


=) 


(ii) 


(111) lim> 


fin £24.9%) - SO) 


6x0 Ox 
— + 6x,t) — F(x, t) 
éx 0 Ox 


F 
Lae + dx,t) — a t) 
x Ox 


6x0 Ox 
Exercise 3 
*F OF @F 0° F 
_ Find 7 Byox’ axdy and ay? for each of the following functions F(x, y). 
(i) F(x, y) = xy 


(ii); F(x, y) =x? + 2x?y + y? 
(iii) F(x,y)=e* 


Exercise 4 


‘Verify that each of the following functions is a solution of the corresponding partial 
differential equation. (f and g are arbitrary differentiable functions of one variable, and a 
and b are given constants.) 


(i) 


eu gu ou 
U = f(x)g(y); axdy ox dy 
(11) U = f(ay — bx); ee 
; ox oy 
= 7 | Gy ay 
(ili) y = f(x — t) + g(x + 2); a ae: 


Exercise 5 


Determine which of the following partial differential equations are homogeneous linear, and 
which are constant-coefficient. 


(i) 


OF 07°F 
x1 + ae = 0 


06°F OF OF OF 


(i) xy->+x-—+y)—->5= 


(iii) x? | 


ax?" dx dy > oy 
OF 4a OF  @F 


cx: py 


OF OF 
(iv) re es ae (where a is a real constant) 


[Solution to Exercises 1-5 on p. 35] 


Summary of Section 1 


1. 


In this séction you have seen examples of partial differential equations and their 
solutions, and you have seen that such solutions involve arbitrary functions, as 
opposed to the arbitrary constants which occur in the solution of ordinary 
differential equations. 


. In this unit we shall restrict our attention to homogeneous, constant-coefficient, 


linear, second-order partial differential equations; i.e. equations of the form 
.6°U 6?U 6*U oU oU 


et aa Ge et 


where a, b, c, d, e and f are constants. 


- Theorem 1: If u and v are two solutions of any homogeneous, linear partial 


differential equation, then for any constants A and B, the function Au + Bv is 


_ also a solution of the equation. 
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2 Properties of coiled springs and elastic 
strings 
2.1 Hooke’s law 


A coiled spring is a device, usually made of steel, which can be altered in length, 
but which returns to its original shape when any external forces are removed. Such 
a spring may be extended, in which case we say that it is in a state of tension, or 
possibly compressed, which we call a state of compression. Our objective in this 
section is to model the behaviour of such springs in some detail, and in order to 
do this we shall need to make certain simplifying assumptions. 


The springs which you have met previously in this course were assumed to be ideal 
springs, that is to. say, they were assumed to obey Hooke’s law and to be light. (By 
‘light’ we mean that the mass of the spring can be ignored because it is small 
relative to the other quantities under discussion.) You will recall that Hooke’s law 
states that an ideal spring when uniformly extended will exert a force F, given by 


where K is the stiffness of the spring, /p is its natural length and / is its extended 


k 
length. We shall find it convenient to replace K by 2 and the constant k is then 
0 


known as the modulus of elasticity of the spring, or more briefly its modulus. We 
then have 


F = ta = lo). 
lo 


The springs under consideration in this section are no longer assumed to be light, 
but I shall make the assumption that they continue to obey Hooke’s law. (This is 
an assumption which can be validated only by experiment.) In future in this unit 
all references to ‘springs’ are intended to mean heavy springs which obey Hooke’s 
law. 


Later in this section I propose to model the behaviour of a spring which 1s 
originally stretched horizontally between two fixed points. The spring is then 
disturbed slightly from its rest position, and our task will be to predict its 
subsequent motion. The essential point to note at this stage is that our starting 
point could well be a spring which has already been extended, but before I 
continue with this line of thought I must define the term tension in the spring. 


We consider a stretched spring with ends A and B which may, or may not, be 
stationary, and consider an arbitrary point P along its length (see Figure 1). 

The section of spring to the right-of P exerts a force towards B, of magnitude T 
say, on the left-hand section; and from Newton’s third law we deduce that the left- 
hand section applies an equal and opposite force on the right-hand section. The 
magnitude of this force is called the tension in the spring at P. The value of T will 
in general depend upon the choice of the point P, but we shall assume throughout 
our discussion that the spring is always stretched rather than compressed, so that 
T is always positive. 


My first task is to examine the relationship between the tension in the spring and 
its physical shape. 


The springs which I shall consider will be long compared to their width, so that in 
future diagrams I shall represent them by straight lines. 


Suppose that the spring is now stretched between two fixed points A and Ba 
distance L apart, and that the spring is stationary. This situation will in future be 
known as the reference state of the spring. 


While in its reference state the tension at all points of the spring is the same, Tp 
say, for each small element of the spring must experience equal and opposite forces 
at its ends, and therefore the tension must be the same at both ends of the 
element. (If this were not so, the element would begin to move in the direction of 
the greater force.) 


Sommmemialtll 
A AWWW IV iB 
P 

AWWW AWA A 8B 

1 
Figure 1. The forces acting at 
P on the left-hand and right- 
hand sections of a stretched 
spring. 


Ty Ty 
| “i — 
A——_—_»—__»—______—— 8 
'¢ Q 


Figure 2. The forces at P and 
Q must be equal. 
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We now choose an arbitrary point P on the spring (and mark it say with a white © 
dot), and our objective will be to describe how P moves when the spring is 
disturbed. Suppose initially that while the spring is in its reference state the point 
P is at a distance x from A. Suppose also that Q is a point close to P, ata 
distance x + dx from A while the spring is in its reference state. | 


The element of the spring between P and Q of length 6x had an original 
unstretched length of da, say, so that applying Hooke’s law to this small section of 
spring we obtain 


rd 


Ty =k 
’ da 


We can solve this equation for 6a and thus find an expression for the unstretched 
length of the element PQ in the form 


da = | ) x (1) 


To +k 
‘Let us now consider what happens if the spring is set in motion by a small 
disturbance along its length. Suppose that at time t the point P on the spring has 
moved so that its distance from A is s, and further suppose that Q is a distance 

s + ds from A. Using Hooke’s law again we have 


: 
5a (2) 
since the length of the element PQ is now 6s, and T is the tension at P. (This 
equation is approximate since we assume that the spring is uniformly extended 
between P and Q in order to apply Hooke’s law.) 


Now, eliminating 6a between Equations (1) and (2), we obtain 
en 
few +h Hk 
Ox 


so that in the limit as 6x > 0 we eliminate the local effects of non-uniform 
stretching of the spring and obtain an exact result 


T= (Ty +k) —k. (3) 


(For the time being I am regarding t as fixed, so that it is unnecessary to use a 
partial derivative in Equation (3).) 
ds 


In Equation (3) you should notice that sf is a quantity which depends upon the 


way in which the spring is distorted, so that this equation relates the geometric 
shape of the spring to the tension at any point of the spring. Also I have not 
actually used the condition that A and B are fixed, so that Equation (3) would 
also apply in the case where A and B are free to move. The above equation is used 
in the following example, which is intended to get you accustomed to the idea of 
tension varying along the length of a spring. 


Example 1 


A stationary uniform spring of mass My, modulus k and unstretched length / 
hangs vertically and supports a weight of mass M. Find the length of the stretched 
spring. 


Solution 


We are told that the spring is uniform, which indicates that its mass per unit length 


M 
when in its unstretched state is constant and equal to oe (When stretched, 
0 


however, the mass per unit length is no longer constant because the stretching is 
not performed evenly along its length.) The solution now falls into two parts: an 
application of Newton’s second law, then an application of Hooke’s law. 


reference state 
A Pu B 
J 
| xX | Ox \ | 


| eae 
eae | 


Y Gee. PdsO B 
a nes 
disturbed state 
Figure 3. The reference state 


and a disturbed state of a 
stretched spring. 
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(i) Apply Newton’s second law to a small element of the spring 


We choose the variables as shown in Figure 4, and consider the forces acting on a 


small element os of the stretched spring. 


The left-hand side of the figure shows the spring in its reference state (with gravity 


‘switched off’) which corresponds to the case in which 7p = 0 and the spring is 
unstretched. This element was originally of length 6x, and its mass is therefore 
Mo 

— 

Suppose now that the tension in the stretched spring is Tp at the point P and Tg 
at the point Q. The element PQ experiences three forces (see Figure 5): 


(1) its weight downwards, 
(ii) the tension at P upwards, 
(iii) the tension at Q downwards, 


giving a net downwards force of 


Mog 


lo 


To — Tp - Ox. 


Since we are assuming that the stretched spring is stationary, the acceleration of 
the element PQ is zero, and from Newton’s second law we have 


force = mass x acceleration, 


so that 


M 
08 ox = Q, 
lo 


(oe 


We now write the difference in tension between Q and P as OT, ie. 
It follows that 


Of Mog 
ae lo é 


and in the limit, 


* eee Se 


Now integrating with respect to x, 


M 
ae ak, 
0 


where C is a constant. 7(x) is thus the tension at the point P which was at a 
distance x from A when the spring was in its reference state. 


In order to determine C we consider the forces acting on the mass M: 
Mogl 
(i) a force upwards of Tg = T(lo) = ee +C=C-— Mog, 
0 
(ii) a force downwards of Mg; 
so that the net force downwards is 
Meg + Mog ee c. 
and this is equal to zero since the mass M is stationary. It follows that 
C=(M + Mo)g. 


(You should notice in passing that we could also have found the value of C by 


examining the tension at the point A and equating this to the total weight hanging 


from A.) 


Here, unlike our previous 
situation, the ‘unstretched 
state’ is our ‘reference state’. 


oe 
unstretched 


spring 
(with gravity 
‘switched off’) stretched 
spring 
(with gravity 


‘switched on) 


Figure 4. The spring in its 
unstretched and stretched 
states. 


Tp 
P 
ds 
QV 
Mos ox To 


li 


Figure 5. The forces acting on 
PQ. 


Tp (= T(h)) 
M | 


Mg 


Figure 6. The forces acting 
on the mass M. 
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Equation (5) now becomes 
Mogx 
T(x) = -— = + (M + Mo)g. (6) 
1) 


Thus Newton’s second law tells us how the tension T(x) depends on x. 
(ii) Apply Hooke’s law 


From Equation (3) on page 10 we have, since Ty = 0 in the reference state we have 
chosen, 3 


dx 


Our problem is to find the length of the stretched spring, which is s(Io) = / say, 
and this would be solved if we could find the function s(x). Eliminating the 
unwanted variable T(x) between (6) and (7), we obtain 


T(x) = K(E -1). 0) 


k fe 2 } a 5 + Male: 
xX lo 


and solving for = we have 


ds M M+M 

eee aa 08 se (M + Mo)g < 
dx klo k 

Integrating with respect to x, we obtain 


5 Mogx? (M + Mo)g 
Lite k 


I: 


s= +i)x+c, 


for some constant C,. 
But s = 0 when x = O, and therefore C, = 0, so that 
_ Mogx? | [(M + Mog 
=— 4+ |. I ] x. 8 
cia | rae. ) 
We are asked to find s(/p), so substituting x = /) into (8) we have 
Mogls . [(M + Mo)g 
— J + IY! 
Real RR es 
(e + 4Mo)g 
k 


S([o) = 
+ i) 


When you derive an answer to a problem like this it is often useful to check its 
validity for extreme cases. Thus, for example, we see that if k is very large, 
corresponding to the case of a very stiff spring, then from Equation (8) we obtain 
_ § =~ x, which is just what we should expect. 


2.2 Vibrations of a stretched spring 


We now come to an important point in this unit, namely the point at which we 
define the problem which we shall eventually solve. 


We consider a heavy spring, at first in its reference state, stretched horizontally 
and tightly between two fixed points A and B at a distance L apart, and then 
disturbed by small longitudinal forces which impart a small motion to the spring. 
More specifically we make the following assumptions. 


(i) The motion of any section of the spring is small compared to the spring’s 
| length, takes place entirely along the length of the spring, and the end 
points are fixed. 


(11) All frictional forces, such as that due to air resistance, may be neglected. 


(111) The forces due to the tension in the spring act along its length and obey 
Hooke’s law. 
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(iv) The gravitational force on the spring is small relative to the tension, and 
so may be neglected. 

(v) The diameter of the spring is small compared to its length, and so may be 
neglected. 


Suppose that we choose an arbitrary point P of the spring, whose reference 
position is at a distance x from A, and suppose that the spring is disturbed slightly 
and begins to move. Our objective is to determine the position of the point P at 
any subsequent time. More specifically, suppose that the point P has moved to a ° 
position at a distance s from A at time t, then s = s(x, t) is some function of x and 
t, and our task will-be to determine this function. Suppose also that the point Q 
close to P, whose reference position is at a distance x + 6x from A, has moved to 
a position at a distance s + 6s from A at time t (Figure 7). 


Before the spring began to move it was in its reference state stretched uniformly 
between A and B, and thus its mass per unit length was a constant, m say. (This is 
not equal to the unstretched mass per unit length.) 


In order to model the motion of the spring we need to bring together two 
ingredients, as in Example 1: 


(i) Newton’s second law (which specifies the acceleration if the force is 
known), 
(11) Hooke’s law (which determines the tension). 


We let Tp be the tension at P and Tg be the tension at Q, at time t. 
(i) Apply Newton's second law to a small element of the spring 


The small element of the spring between P and Q was originally of length 6x, so 
that its mass is mdx. This element is acted on by two opposing forces, the tension 
at P and the tension at Q. The centre of mass of this element is somewhere 
between P and Q and, since the element is small, its acceleration is approximately 

67s 

Ot? 
(This corresponds to the acceleration of a particular point of the spring, i.e. a fixed 
value of x, and hence we keep x fixed in the differentiation.) From Newton’s 
second law we have 


(x, t) | to the right. 


force = mass x acceleration, 


so that 
075 
To = Tp ~ mox Or” 
and therefore 
6T . 0*s 
eee: soak 
In the limit as 6x — 0 the centre of mass of the element becomes P, so that 
6T 0*s 
eee 9 
ax ar 0) 


(ii) Apply Hooke’s law 


We now use Equation (3) on page 10 (derived from Hooke’s law): 
ds 
T = (To + kK) —k, 


The derivative in this equation refers to a fixed value of t. Now that we are 
ee Se ; , 
treating t as a variable we must write a to indicate that tf is held constant during 


the differentiation, and we then have 


Q 
Pati. + k) ok, (10) 


reference state 


A PQ. B 
pr eae enema 
l x | dx \ l 
| La | 
| ze | 
Ae FOC iB 
fe 
disturbed state 


Figure 7. The spring in its 
reference state and a disturbed 
state. 


SEE EES ane EON 
Pe oO 


Figure 8. The forces acting on 
a small element of the spring. 
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Differentiating (10) partially with respect to x, we obtain 


OT 
a 
We wish to discover how s depends on x and ft, and therefore we want to eliminate 


the superfluous variable 7. We can do this by equating Equations (9) and (11) to 
obtain . 


| 07s 
(To + Ka (11) 


0*s 6*s 
To tk. = a ; To +k 
The value of a is positive, and it is common practice to put c? = ee to 
produce the partial differential equation 
(13) 


This equation is known as the (one-dimensional) wave equation, and is of crucial 
importance for the rest of this unit. In fact it is this partial differential equation 
that you will learn how to solve. The constant c is called the wave speed. 


The function s = s(x, t) determines the distance of a point P from A at any time [; 
but it is more common to write the equation in terms of U = s — x, so that 
U = U(x,t) determines the distance of P from its starting point (see Figure 9). 


If we substitute s = U + x into Equation (13), we obtain 


‘PU 1 @U 


= 14 
ay ee se 


Exercise 1 
Verify that the substitution of s = U + x into Equation (13) gives Equation (14). 


[Solution on p. 35 | 


The advantage of using U rather than s is that the conditions which follow from 
the assumption that the ends of the spring are fixed are quite simply stated as 


U(0,t) =0 
for all t. 
U(L, tf) = 0 


These conditions are called the boundary conditions for the partial differential 
equation, and will be used extensively in Section 3. 


2.3 Elastic strings 


An elastic string is just what it says it is, a string with elastic properties; but as a 
mathematical term it specifies that the string satisfies Hooke’s law. Springs and 
strings are therefore identical in behaviour except in one important respect; a 
string cannot be in a state of compression. 


We introduce the concept of an elastic string simply to give you another 
illustration of the use of the wave equation, and one which is perhaps a little easier 
to visualize. We suppose that the string is initially stretched tightly between two 
fixed points A and B, and then disturbed slightly in a direction perpendicular to 
the length of the string (see Figure 10). 


If P is a point on the string at a distance x to the right of A, then it can be shown 
that U(x, t), the perpendicular displacement of P from the line AB, will satisfy the 
wave equation 

6*U a8 6*U 


ax? sc? Ot? 


reference state 


disturbed state 


Figure 9. Derivation of the 
function U(x, fr). 


Figure 10. An elastic string 
vibrating. 
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where c is some constant, dependent on the physical properties of the string and 
its initial tension. (I ask you to accept this on trust since I feel that it would 
take too much of your time to give yet another derivation of the wave equation.) 


In this case the boundary conditions are again 


U(0,t) =0 


Ee: 
sagen all t 


As in the case of a spring, the task is to determine the function U(x, t) from the 
initial behaviour of the string. 


In order to distinguish the two forms of motion, vibrations along the length of the 
string or spring are called longitudinal vibrations and vibrations across a string at 
right angles to its length are called transverse vibrations. 


Summary of Section 2 


1. In this section we modelled the behaviour of a ‘heavy’ spring, and we used 
Hooke’s law to derive the equation 


d 
{7 + = = 


2. By combining this result with Newton’s second law we showed that with certain 
assumptions the longitudinal displacement U of a stretched spring satisfied the 
wave equation 


6*U et 6?U 
ose? oF 
for some constant c, where U satisfies the boundary conditions 


U(0,t) = 0 


for all t. 
ui ve 2 os 


3. We introduced the terms longitudinal and transverse vibrations, and noted that 
transverse vibrations of a string also satisfy the wave equation and the same 
boundary conditions. 


3. The method of separation of variables 


You will recall that I have set myself the task of finding a function which 
determines the position of a vibrating spring at any time t¢. In other words, I have 
to find a specific solution of the wave equation 

ioe | ee Si siae 8] 


ie eS or 
satisfying the boundary conditions 
U(O.t) = 16, b= 0 for all t. 


It is clear that this equation as it stands does not contain sufficient information, 
for at no stage in the previous section have I introduced the initial position and 
velocity of the spring. However, as we shall see, if the initial position and velocity 
are specified we can find a definite solution to the problem. 


Before we get deeply involved in the details of the methods which we use to solve 
partial differential equations, I should like to say a little about the general strategy 
that we adopt. The point is that in general it is really rather difficult to find all the 
solutions of a partial differential equation. If, however, we restrict our search for 
solutions to those of a particularly simple form, then we may be able to reduce the 
problem from a difficult partial differential equation to a number of much easier 
ordinary differential equations. Such a process is the method which is known as 
separation of variables, and | shall explain the details of it shortly. 
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Let us suppose for a moment that we do have a method which will enable us to 
find some of the solutions of the wave equation 

6?U 4 6°U 

oe gt oF! 
say U,, U2,...,U,; then since this partial differential equation is homogeneous 
linear, we know that any linear combination 


», AU; 
i=1 


(where A,, A2,...,A, are constants) is also a solution, from Theorem 1 on page 7. 
The method of separation of variables does indeed provide us with such a set of 
solutions, and by taking linear combinations we can obtain many more. But this is 
an idea which we shall develop later; at the moment our concern is to find some 
solutions of the wave equation. 


3.1 Separation of variables 

Our immediate plan is to find some solutions of the wave equation 
6*U 4 6*U 
Oe oF 

satisfying the boundary conditions 
U(0,t) =0 
L(L, =O 


We now take ‘a leap in the dark’, guided only by the fact that other 
mathematicians, wiser than ourselves, have seen the final outcome, and attempt to 
find some solutions of Equation (1) by assuming that they take the form 


(1) 


for all tf. 


U(x, t) = X(x)T(t). (2) 
In other words. U is the product of two functions X and T. where X is a function Do not confuse this use of T, 
of x alone and T is a function of t alone. (There is of course at this stage no to indicate a function of time, 
: ; with our previous use of this 
guarantee whatever that there are any solutions which take this form. By ide ta ennie (ee tersion in 
constructing them we shall in fact show that they exist.) a spring. 


Differentiating Equation (2) partially with respect to x, we obtain 


UG) _ x'@)T() 
Ox 


where X’ denotes the derived function of X. (Remember that T, being a function of 
t only, can be regarded as a constant for this calculation.) Differentiating partially 
with respect to x again, and abbreviating the notation, 


6°U 


al tee yy’ T 

Ox? (3) 
Similarly, differentiating the original Equation (2) twice with respect to t, we obtain 

0*U 

md F 

sy =X1 (4) 

Substituting Equations (3) and (4) into (1), we get 
” 1 ” 


Now dividing both sides by XT to get all the functions of t on one side and 
functions of x on the other, we obtain 

f Ss = Dea 

Cf =X 
(assuming that neither X nor T is zero). This then is a necessary condition that 
U = XT should be a solution of the original equation. 
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The next step caused me considerable difficulty as a student, and I needed to think 
about it for quite a while. We let 
= F aus = 7 
ST ¥ 


” 


then since yp = (a function of t only), it must be independent of x; but since 


c7T 
ut 


H= > (a function of x only), it must be independent of rf. But if uw is independent 


of x and t, it can only be a constant. 


[I hope that I have made the point more clearly than it was made to me, but, if 
4) = X"(x) 


you still find it difficult, let u(x,t) = oT eae Then 
gx, = Fir = (0, t) (since each expression is independent of x) 
“ 
and 
x “(x) er ee 
u(0, t) = X (x) = 1 O-4)) (since each expression is independent of f). 
Hence 


u(x, t) = u(0,0) for all x and f, 
and therefore y is constant. | 


It follows that any solution of the original partial differential equation which is of 
the required form must satisfy the two ordinary differential equations 


T"=cyT and X"=pX. (5) 


We know that yi is a real (rather than a complex) constant (because T and 7” are 
real), and therefore there are three cases to consider: 


gS 0, =O ane aa 


Before I discuss the three cases in general, try the following exercises, the first to 
prepare you for the general argument and the second because the result will be 
useful later. 


Exercise 1 
Solve the equations 


7 Secul and X" = px 
in each of the following cases: 
tes ee 
(i) c= 22. p= 1 
tii) c= — 2, p= 1b 


Exercise 2 
Ifh £0, and Ae” + Be~™ = 0 for all values of t, show that A = B = 0. 


[Solutions to Exercises 1 and 2 on p. 36 | 


Remember that we are trying to solve the partial differential equation 

6*U 4 6°U 

pe oe Or 
with U(0, t) = U(L, t) = 0, and we have seen that a solution of the form 
U(x, t) = X(x)T(t) must satisfy T” = c?u T and X” = » X. Now we use the 
boundary conditions to filter out the solutions which are of no interest in the three 
cases p > 0, uw = 0 and uw < 0. 


The process is rather like panning for gold — we need to filter out quite a lot of 
useless material, until finally we are left with a nugget at the bottom of the bowl. 
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Case I: up > 0 

For convenience we can put p = h? where h is a real constant; then the equations 
TT’ acetone X" = h’Xx 

have solutions 
T=Ae™+Be™ and... 4% = Ce% ade 

where A, B, C and D are constants. 


Now if either T(t) = 0 for all t or X(x) = 0 for all x, then U(x, t) = 0 for all x, t, 
and although this is in fact a solution of the partial differential equation, it is a 
rather boring solution since it corresponds to the spring remaining stationary. The 
solution U(x,t) = 0 is called the trivial solution. We are actually looking for non- 
trivial solutions which satisfy the boundary conditions 

U0; t) = UfLa) =0 for all t. 
These boundary conditions refer to particular values of x (namely 0 and L) and all 
values of t, and assert that 

(Ce°. + Dey Fit) =0 

(Co + De ™ Tit) = 0 
As we are searching for non-trivial solutions, we may assume that T(t) 4 0 for 
some value of t, and this allows us to conclude that 

Cc +P =(° 

Ce + De =v 
The Unique solution to this pair of equations is C = D = 0, which leads to the 
trivial solution after all, in spite of all our efforts to avoid it! 


| for all t. 


Case II: un = 0 
In this case the equations for T and X become 


i «= and x ee ©: 


As in Case I, all that really matters about T is our assumption that T(r) 4 0 for 
some t. The solution of X” = 0 1s 


A = Ox, 
and if we again apply the boundary conditions we obtain 


_ DIT(t)=0 
(CL + D)T(t) = 0 


Our assumption that T(t) 4 0 for some t implies that 
Def 
CL+D=0 


; for all f. 


so that C = D = 0, and again the only permissible solution is the trivial solution. 


Case III: uw < 0 
For convenience we put up = —h’; then 
X"=-Wx and Fos er. (6) 


As in part (ii) of Exercise 1 we see that the solution of the first -of these equations 
is 


X = Ccoshx + Dsinhx. 
We now apply the boundary conditions (again simply assuming that T is not the 
zero function) to obtain 


CT(t)=90 


(C coshL + DsinhL) T(t) = for all 1, 
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and then choose a value of t for which T(t) 4 0. We deduce that C = 0 and 
therefore 


DsinhL = 0. 


We may assume that D # 0 (for otherwise we are once again left with the trivial 
solution), and therefore 3 


sinhL = 0. 


This is the ‘gold nugget’ for which we have been-searching. If sin hL = 0, then it 
follows that 


hL=ae 18. ae 
z 

for some integer n. 
Solving the second equation of (6), we obtain 

T = Acoscht + Bsin cht, 
which gives us the solution 
U(x, t) = X(x)T(t) 

= Dsinhx(A cos cht + Bsin cht) 


nt 
where h = = for some integer n. 


Substituting this expression for h and absorbing the arbitrary constant D into the 
arbitrary constants A and B, we finally have a solution 


ex cntt cntt 
U(x, t) = sin — [4 


3 er a + Bsm—— 


which satisfies the boundary conditions and which is certainly not the trivial 
solution. More than that, we have discovered an infinite number of distinct 
solutions corresponding to the choices of the integer n: 


aaa cmt _ cmt 
U,(x,t) = eet [a weet a + B, sin) 


eX 2crt <- 2ent 
bix.t) = oa I. ae Se + B, ace 


etc., and the general function of this type would be 


Te sea nx oe nent ha. ncnt 
— ——— <= he 
al Xs sin 7 , COS 7 nS 7 


where n is an integer, and the As and Bs are arbitrary constants. 


These solutions resemble the normal modes of a lumped-parameter vibrating Unit 24 
system, and what we shall do next corresponds to taking linear combinations of 

normal mode solutions to get the general solution. Although we have indeed found 

rather a lot of solutions of the wave equation, we have by no means found them 

all, and it will need some further work to increase our store of solutions. 


Exercise 3 
Use the method of separation of variables to show that if 


A(x, t) = X(x)T(t) 
is a solution of the partial differential equation 
ee | eT. 
= ¢€ Soa 
ax? ot 


then 
Aa ak ond T= pt 


for some constant p. 
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- (This partial differential equation may be used to model the flow of heat along a thin metal 
rod, in which case 6(x,t) represents the temperature at time t at a distance x from one end 
of the rod.) 


[Solution on p. 36] 


3.2 Initial conditions 


In the last subsection we discovered an infinite number of solutions of the wave 
equation (each satisfying the appropriate boundary conditions): 


nict . nnct 
—— + B, sin —— 


nTx 
U(x, t) = sin — |A, 
(x, t) = sin L = L L 


where = 1, 2 32. 
The current state of our investigation is summarized in the set diagram, where: 


SF = {all solutions of the wave equation}, 


% = {all solutions satisfying the boundary conditions U(0,t) = U(L,t) = 0}, 


FP = {all solutions of the product form U(x, t) = X(x)T(t)}. 


We have discovered that Pm &, the intersection of the sets Y and Z%, consists of 
solutions {U,(x, t)$; but our objective was to discover all the elements of the set Z 
(shaded in Figure 1). 


We now use Theorem 1 (page 7), which tells us that we may add two solutions of 
the wave equation and obtain a third solution, but there is still the question of the 
boundary conditions. In fact if we add two solutions which satisfy the conditions 


U(0,t) = U(L,t) = 0 for all t, 


then we obtain a third solution with this property. For this reason these particular 
boundary conditions are called homogeneous. (We shall see an example of non- 
homogeneous boundary conditions in Section 4.) 


In fact any function of the form 


N 


U(x,t) = Y U,(x, 0) 


n=1 
for some integer N, is an element of the set Z. 


We are now ready to introduce the final factor into our investigation: 


The initial conditions 
In order to specify the motion of the spring, we need to know where it is and how 
it is moving at the start of the motion, say at time t = 0. To know where it is we 
need the function U(x,0), which tells us the displacement of each point of the 

Se ieee ©) 
spring when t = 0. To know how it is moving we need the function 3p 09), 
which tells us the velocity of each point of the spring when t = 0. 


If we are given U(x,0) and —(x,0) then this should be sufficient information for 


us to pick out the function U(x, t) from all the other functions in the set BZ, and 
thus to determine the motion of the spring. Let us now see how this works in 
practice. 


Example 1 
A uniform heavy spring stretched horizontally between the points x = 0 and x =z 
is given an initial displacement 

U(x,0) = sin x 


along the length of the spring and released from rest. Find its subsequent 
displacement as a function of x and t, if the wave speed (i.e. the constant 
appearing in the wave equation) is c. 


Figure 1. Various sets of 
solutions of the wave 
equation. 


The connection between the 
two senses of the word 
‘homogeneous’ (as applied to 
partial differential equations 
and as applied to boundary 
conditions) is that they both 
imply the property that ‘the 
sum of two solutions is 
another solution’. 
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Solution 

We can apply the wave equation 
6*U 4 6*U 
at ir 


for some constant c, and we try a solution of the form 


Utx,t}= y UO Ax) 


N 
= ) sinnx(A,cosnct + B, sin nct) (7) 
n=1 
for some integer N, where A, and B, are constants which have yet to be Note the simplification that 
determined. . results in this example from 
choosing the length L to be 
We are told that equal to z. 
U(x,0) = sin x, (8) 
and since the spring starts from rest we have 
oU 
——(x,0) = 0. ? 
5 (0) (9) 
Putting t = 0 in (7), and using (8), we obtain 
sinx = A, sinx + A,sin2x + A3;sin3x + --: + Aysin Nx. (10) 
Then differentiating (7) partially with respect to t, putting t = 0 and then using (9), 
we obtain 
0 = c(B, sinx + 2B, sin2x + 3B3sin3x + -:- + NBysin Nx). (11) 


Our task is now to choose the As and Bs in order to satisfy Equations (10) and 
(11), and I have chosen this example deliberately so that it is possible to see these 
values at once. All we have to do is to choose 


A, = 1, A, = 0 for n> 1, to satisfy Equation (10), 
and 

B, = 0 for all n, to satisfy Equation (11). 
Equation (7) now reduces to 

U(x, t) = sin x cos ct 
and this is the solution to our problem. 


At different instants of time t, the function U(x, t) tells us the displacement of the 
point on the spring whose reference position is at a distance x from A. 


If we choose a particular point P on the spring, so that x is regarded as constant, 
then P moves backward and forward under the influence of the factor cos ct. The 
amplitude of the vibration is determined by the factor sin x, so that the movement 
is, as we should expect, most violent at the centre of the spring. 


Figure 2 (overleaf) shows the positions of various points of the springs at various 
times; for example P,, P,,...,P; represent various positions of the moving point P 
at times. t,t; ... ; fs. 
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reference state 


initial displacement 


Figure 2. The movement of a longitudinally vibrating spring. 


It can be seen from this example that our task is finally reduced to finding the 
coefficients in two expressions similar to Fourier sine series, and thus at last we see 
the link between this unit and the preceding one. Using the techniques developed 
in Unit 31, we should be able to solve our problem in general. 


3.3 Finding the coefficients in the Fourier series 


Suppose that ¢ and yw are given functions, each defined between 0 and L, with the 
_ properties that 


(0) = P(L) = WO) = W(L) = 9, 


and we wish to find a solution of the partial differential equation 


CN CU 
ox ee 
on the interval 0 < x < L, which satisfies the boundary conditions 
U(0,t) = U(L,t) =0 for all t 
and the initial conditions 
U(x,0) = d(x) TOs x < LL. (12) 
—(x,0) = (x) f0<2<2 7. (13) 


Using the results of Subsections 3.1 and 3.2, we can attempt to find a solution in 
the form of an infinite series 


oC 


Uist) = > Uae 


, n=l 
SO nict : 
= ‘3 sin A, cos——~ + ae Sooners (14) 


n=1 
From Equation (12), and putting t = 0 in the above series, we obtain 


o(x) = U(x, 0) = y Assit (15) 


n=1 
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From Equation (13), and differentiating the series in Equation (14) partially with 
respect to t, we have 


oU Say NTC nuct Nnmc nict 
——(x,t) = in —— | - —A, sin + ——B, 
(x, t) 2 sin = sin te L 
and therefore 
w(x) = “(x,0) = ee (16) 
n=1 L L 


We now use the techniques of Unit 31, Section 4, to determine the constants A, 
and B,,, and to do this we use the functions @ and wW (which are non-periodic) to 
construct two new functions, ® and ‘V, which (like the sine function) are odd 
functions and periodic. We let 


D(x) = o(x) i Oe xk 
=x) abe ee 


and 
V(x) = wx) if fe 


< 
—w(—-x) if -L<x< 


and further require that both functions ® and PV are periodic with period 2L (see 
Figure 3). 


v= G(x) or Wx) 


Figure 3. Obtaining an odd function of period 2L from a function defined between 0 and L. 


We now use the formula for A, given in Subsection 4.1 of Unit 31 , allowing for the 
fact that the period of our function is 2L, not L. We obtain from Equation (15): 


Amz D(x) sin—— dx 


=F, 
=i (fo x) sin —~ dx = | (x)sin ™* dy] 
0 L 


Se a _ AMX 
or -(f o(x)sin —— dx + I #(—»)sin"™ dy] 


Sas _ ATX : _ AT 
rs if (x) sin—— dx + i} g(t) sin as 


(putting : = —x in the second integral). 


Thus 


a _ AMX 
A, = al (x)sin—— dx. 


Similarly we obtain from Equation (16): 


—- 


oa B=i W(x)sin— dx, 
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2 
B, = — 
ATCC 


L 
{ W(x) sin a dx. 
0 L 


To summarize, the general procedure for solving the wave equation with 
homogeneous boundary conditions is given in the box below. 


Procedure 3.3 
Given the wave equation 
ou 490 


a <x< 
2 


subject to the boundary conditions 
U(0,1) = UG, =0 for all t¢ 


and the initial conditions 


— (x,0) = W(x) foes L. 
1. Use the method of separation of variables to obtain solutions 


nix nict nict 
U t) = sin —— |A, cos —— + B, sin——]}. 
Re ae 3 | sin (4, 005" + B, SI L 


2. Use the facts that the equation is homogeneous linear and the 
boundary conditions are homogeneous, to construct a general 
solution 


U(x,t)= > sin s cos —— + B, sin ——}. 


= nlx 7 nict —pepeet 
n=1 ‘ L L 


3. Use the procedure in Subsection 4.1 of Unit 31 to obtain 
formulas for the As and Bs as follows: 


L 
A= al (x) sin dx 


ae nmx 
= aise eae aD ESE ay d e 
B,, [ W(x) sin 7 x 


If the functions ¢@ and w are sufficiently simple, we may attempt to find formulas 
for A, and B, in terms of n; otherwise we may compute say the first twenty 
coefficients, perhaps using Simpson’s rule and a modern computer to calculate the 
integrals. In either event our problem is essentially solved. 


Let us try this technique on a particular problem. 


Example 2 
Solve the partial differential equation 
eu _10u 
Ox ee oF 
on the interval 0 < x < z subject to the boundary conditions 
U(0,t) = U@=0 for all t, 
and the initial conditions 
U(x,0) = x(x — x) rl <= x= 4, 


(x40) =0 fOcx<nz. 
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Solution 
We attempt to find a series solution 


U(x,t) = )° sinnx (A, cos nct + B, sinnct). 


n=1 


Applying the formulas in the above box for A, and B,, we obtain 
a T 
An = = | x(x — x)sin nx dx. 
1 Jo 
Using the method of integration by parts, we can show that 


4 
Tl 
{ x sin nx dx = ——cosnn 
0 n 
and 


: 1 3 
| x? sinnx dx = ——cosnn + =; (cosmn — 1). 
0 n n 


Thus 


8 ee 1 2 
A, = | — cos +-——cos nn — —, (cosnx — vf 
T n n n 


4 
= —,;(1 — cosnz) 
n 


4 n 
Set 1h) 
Since the initial conditions give w(x) as the zero function, we have 


5. = for all n, 


25 


and thus we have found a solution to the problem in the form of an infinite series: 


ie. 6) 


4 1 
U(x, t) = = > 3 ll — (—1)")sin nx cos net. 
n=1 


This expression can be simplified, by noticing that 1 — (—1)" is equal to 0 for 


even n and to 2 for odd n, so that we have to sum only over odd n. Thus, putting 


n = 2m — 1, we finally obtain 


8 00 
U(x,t) =— ——— sin (2m — 1)x cos(2m — 1)ct. 
(0) = Gg py sin 2m — 1x cos (2m — 1) 
Exercise 4 
Solve the partial differential equation 
67U ce. 6?U 


6x? c? Or? 
with the boundary conditions 

U(0, t) = U(x, t) = 0 for all t 
and initial conditions 


U(x, 0) = o(x) fhe x<ee 


oU 
— (x,0)=0 fO<x<nz 
ot 
where (with d a given positive constant) 
(x) = if0<x<- 
T Z 


= 2d 


'-*] tte 
1 2 


Hint: the solution to Problem 3(i) in Unit 31, Section 5 shows that for any positive /, 


oe Sd, 
7 |, dtx)sin™ ax = Soin, 


Just as in Example 1, the 
assumption that L= zx 
simplifies the formula. 


Figure 4. The graph of (x). 


26 


where 


2dx 


l 
oma}, : 
2a{ -*] if 


[Solution on p. 36] 


Summary of Section 3 
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In this section we have used the method of separation of variables to obtain a 
solution of the wave equation subject to certain boundary conditions and initial 


conditions. 


The method is summarized in the following flow chart. 
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The problem 


Boundary conditions Wave equation Initial conditions 


U(0.1) = U(L.t) = 0 
U(x,0) = d(x) 


for all ¢. : 
0 

= (x,0) = U(x) 
or 


Look for solutions of the form 


O(x.t) = X(x) T(t). 


Separation of variables gives 
PF =7-a7 ond =X 


for some constant p. 


Yes 


p=? 


U(x.t) = (Cx + D)T(t) 
where T (ft) remains 
to be found 


T=Acoscht + Bsin cht 


X = Ccoshx + D sin hx 


U(x,t) = (Ce™ + Dew") T(t) 


gig where 7 (ft) remains to be found 
The boundary conditions imply 


C=(OQandDsinhLl =(0, 
Only the trivial 
solution obeys 
the boundary 

conditions in 
these cases. 


giving non-trivial solutions if 


h =nn/L for some positive integer n 


There are solutions 


7 


- Aa 
U,,(x.t) = sin —— | A,, cos 


L L + B,, sin 4 


na7ct - et | 


The general solution with these boundary conditions is 


uy ee ao , niet B <i nnct 
= > i n COS — + n SiN—>— | 


n=1 


Use the fact 
that we can 
add solutions. 


2 
i 


3 iS 


yx) sin —— dx 
0 L 


s 
(x) sin — dx 
0) L 


ATC 
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4 Transverse vibrations of a stretched 
string (Television Section) 


The following section may be read either before or after the television programme 
‘Waves — Partial differential equations’. 


Suppose that an elastic string is stretched tightly between two fixed points A and 
B, and disturbed slightly in a direction perpendicular to the string; then the string 


will begin to vibrate. If U(x,t) denotes the displacement at time t of a point P ee 
which was originally at a distance x from A, then U satisfies the wave equation 

6?U =, 6?U 

Oe eee 
for some constant c. 
The methods which we have developed for solving this equation can be applied 
with no change at all, and it is only in the interpretation of the solution that the 
difference from the case of longitudinal vibrations occurs. 
Example 1 
Suppose that the points A and B are a distance L apart, that U satisfies the wave 
equation, and that the initial conditions are 

U(x,0) = sin if0<x<L, 

oU 

—(x,0) = 0 ris x= L. 

ot 
Then using the same procedure as in Example 1 of Subsection 3.2 (page 20), we i 
obtain P, 

aX et 
U(x, t) = sin— cos —. 
(x, t) L L 
Figure 1 shows the position of the string at various instants of time, where 
P,,P2,...,P.¢ represent the positions of the moving point P at times t,,t»,...,tg. 
Essentially then, this problem of the transverse vibrations of a stretched string is 
completely solved. 
a 

New boundary conditions x : s 


Let us now consider the effect of changing the assumptions in the problem, and S 
instead of fixing B we shall assume that some external influence causes it to Figure 1. The positions of a 
° : : ; : ‘ ‘ ; vibrating string at six instants 
vibrate at right angles to the string in such a way that its displacement at time t is of time. 
Rsin at. 


The function U still satisfies the wave equation, but we have altered the boundary 
conditions to 


U(0,t) =0 
for all t. 


U(L, t) = Rsinowt 


This is an example of non-homogeneous boundary conditions, for if we add two 
solutions of the wave equation which satisfy these conditions then we do not 
obtain another solution with this property. You should work through the next 
exercise to discover how this difficulty may be overcome. 


Exercise 1 

Assume that U(x,t) = X(x)T(t) is a solution of the wave equation, so that 
T" =c*pT 
xX” = ux 


for some constant p. 
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(i) Use the boundary condition 
U(L, t) = Rsinot 
to show that 


Rsin wt 
1i)}= 
(ii) Substitute the above expression for T(t) into the equation T” = c?yT and show that 
z 
a 
ace SE. Te 


Cc 


— 


(i111) Solve the equation 


Ges = 

z 

with w= ——-. 
Cc 


(iv) Use the condition U(0,t) = 0 to simplify the solution obtained in part (111). 


(v) Write down a solution of the wave equation which satisfies the above (non- 
homogeneous) boundary conditions. 


[Solution on p. 36] 


At this stage you may well ask how we might find the many other solutions of the 
wave equation which also satisfy these boundary conditions. In fact we use a 
rather nice trick, and suppose that U (x, t) is the solution which we have found in 
Exercise 1, and U(x,t) is some other solution satisfying these boundary conditions. 
Then consider the function 


V(x, t) = U(x, t) — U(x, t). 


The essential point to notice is that Up and U satisfy the same condition at x = L, 
namely 


U{L, 1) = U,{(L,4) = Rein ot for all t, 
and therefore 
V(L, t) = U(L,t) — UQ(L, t) = 90 for all t. 


Thus V(x, t) is just a solution of precisely the problem which we solved in 
Section 3, i.e. it satisfies the wave equation and the boundary conditions 


V¥O.0 = ViL,t) = 0 for all t. 

Thus we can find V(x, t) as a series and then write 
U(x, ft) = 8 oix, t) + Vex, 0) 

to solve our present problem. 


An example should make the method clearer, and for simplicity I shall choose 
t= 7. 


Example 2 

Solve the wave equation 
6*U Ss 6*U 
ax Cer 


with boundary conditions 
U(0,t) =0 for all t 
U(x, t) = Rsin wt for all t 
subject to the initial conditions 
U(x, 0} =0 for 0:< x <= 4, 


(x0) =0 forO <x <7. 
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Solution 
We begin by letting 


ee ee 
Rsin — sin wt 


Z 
U(x, t) = a 
sin —— 
Cc 


which we know from Exercise 1 satisfies the wave equation and the boundary 
conditions but not the given initial conditions. 
Then we put 

V(x, t) = U(x, t) — Uo(x, t). 


Putting t = 0, and using the notation of Subsection 3.3, we have 


V(x,0) = $(x) 
where 
p(x) sts U(x, 0) si U o(x, 0) 
= forO <x <u. 
Also, 
OV 
(50) = Wer) 
where 
oU el OF 
W(x) pes ew 5 at (x, 0) 
Rw sin — 
ais 
= eR 
sin —— 
£ 


Using the box on page 24, we have 


co 


V(x,t) = > sinnx(A,cosnct + B, sin nct) 


n=1 
where 
A, =0 for all n (as ¢ is the zero function) 
and 
2 TT 
B, = — | w(x)sinnx dx 
NIMC Jo 
—2Rw + pe 
———————— | $i ——- 6H AX dx 
. WD 0 
nic sin —— 
40) Tt 
sin !|— —n]x sin |—+n|x 
— Rw C 
Ww Ww 
nic sin —— ——7 == 
Cc Cc 0 


WM 7 
2 + nls 
ra 
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— Rw wl. [@ . lw _ lo . |@ 2 
= ————— | —|sin|— — n]x — sin|— + n]|x] + nI{sin|j— — n}]x + Ssinij— + n|x 
: Cc e c c 
nic sin —— 5 
—}| —n 
6 0 


wx 
et cos — sin nx + 2n sin —— Cos nx 
c 


os. 
2 _ on F ) -" 
nc sin — 
C 0 
oss the trigonometric formula for sin (a + #) on page 8 of the Handbook, letting 
= — B= +n) 
—2Rac 


_ on 
= —__________ |sin — cos nn 
_ on 
n(w? — n*c*) sin — 
C 


= stom ea 
ae — act)” 
Finally we have for our solution for V, 


Roc oe) iy" 
+g Asin nx sin nct 
feel ee n? 


VX, t) = 
and our solution for U(x, t) is 
Rsin ab sin wt 
—sin w 
gee 


sin —— 
Cc 


oo 


oe a nx sin nct. 
_- 


This then is the solution of our problem (which is referred to in the television 
programme), but you should notice that it is certainly not valid if sin —— = 0, in 
C 


other words if @ is some multiple of c. In fact if @ is close to some multiple of c 
then the function 


SX 
Rsin —sin owt 
U (x,t) = 
sin —— 
Cc 


will be very large in modulus for some values of t. One of the terms in the series for 
V(x, t) will also be very large in modulus, and near t = 0 these two large terms will 
effectively cancel out, but as we progress in time, these large terms will dominate, 
and we should expect some very wild oscillations of the string. 


(There is a further difficulty which might have occurred to you. The point for 
which x = x (attached to the vibrator) is at a distance R sin wt from its rest 
position at time t. Its velocity is therefore Rwcos at at time t, and hence when 

t = 0 it has a velocity Rw perpendicular to the string. But all of the adjacent 
points of the string are assumed to have zero velocity. There is no mathematical 
difficulty here, but there is certainly a physical one. To overcome this problem we 
might assume that there is a very small link, between the vibrator and the string, 
which is allowed to move, and whose effect we ignore.) 


The strategy which we have adopted to solve this problem and the physical 
interpretation of the solution are discussed in the television programme which 
accompanies this unit, ‘Waves—Partial differential equations’. 


Summary of Section 4 


The flow chart overleaf is intended to provide a summary of this section. 
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The problem 


Initial conditions 
U(x.0) = 0 


Boundary conditions Wave equation 
PU Aw 


- , 
OX 


ee yt- AU 
Ee — (4,0) = 0 
ol 


INHOMOGENEOUS 
BOUNDARY 
CONDITION 
U(L.t) = Rsin wt 


Look for a solution of the foftm 
UO (x.t) = X(x) T(t). 


Separation of variables gives T” = c- wT 
and X” = wX (constant). 


R sin wt 


T(t) = XL) 


Bring together 
the expressions for 
X(x) and T(t): 

D cancels out. 


= 
Rsin — sinwt 


Up(x.t) = 


The general solution with these boundary conditions is 
U(x,t) = Uv (x.t) + V(x,t) 
where V(x,t) is the general solution of the problem 


with homogeneous boundary conditions. 


Use the flow chart on page 27 
to find V(x,f). 


The general solution with these boundary conditions is 


. WX , 
R sin — sin wt ae 


c _ nm ntct = eater 
U(x,t) = eee be sin >= — (4. cos, + B, sin L 
sin a2 n= 
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5 End of unit problems 


In this set of problems you will be led through an investigation of the vibration of 
a uniform heavy spring suspended vertically from a fixed point with a weight of 
mass M attached to its lower end. Later for simplicity we assume that M = 0. 


We use the notation of Example 1 of Section 2 (page 10), and suppose that the 
spring is of unstretched length /), modulus k and mass Mo, and that a mass M is 
attached to its lower end. 


The element between P and Q is presumed to have moved to become an element 
of length ds at time t, and P is then at a distance s from A. 


You may attempt this set of problems in the following form, or if you find it too 
difficult you may work through the sequence of questions given after the problems, 
which take you through the working in easy stages. 


Problem 1 

With the above notation and assumptions, show that 
6s 4 Mog a Mo 65 
Ox? lo Ss lo ét? 


Problem 2 


Find a solution s(x, t) of the above partial differential equation which is independent of time 
(that is s(x, t) = X(x)). 


Problem 3 
Assuming now that M = 0: 


(1) | show that the tension in the spring at B is zero; 


(ii) | use Hooke’s law to show that 
0 
id) = 1 for all t; 
Ox 


(iii) find the general solution of the above partial differential equation in the form of a 
series, subject to the boundary conditions 


s(0,t) = J 
oar | for all t. 
Ox 


(Hint: Consider the difference between the required solution and the solution which is 
independent of time.) 


[Solutions to Problems 1—3 on pp. 36-8 | 


Sequence of questions 


Our objective is to find a partial differential equation involving s(x, t), the distance 
of P from A at time t, and we shall require the tension in the spring at P, which we 
denote by 7. Read again through Example 1 of Subsection 2.1 (page 10) and then 
answer the following questions. 


(a) | What is the mass of the element ds? 
(b) What forces are acting on the element os? 
(c) What is the total downward force acting on the element 6s? 


(d) Write down expressions for: 


(i) the distance of P from A, 
(ii) the velocity of P downwards, 
(iii) the acceleration of P downwards, 


at time t. 


(ec) Use the answer to Question (d) to write down an equation derived from Newton’s 
second law 


force = mass x acceleration 


C2 
involving T and Take the limit of this equation as 6x > 0. 


—_— 


unstretched spring 
(with gravity 
‘switched off’) 


stretched spring 
(with gravity 
‘switched on’) 


Figure 1. A spring in its 
unstretched and stretched 
states. 
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A 


(f) | Use Hooke’s law to write down an expression for T in terms of ro (See Subsection 2.1.) 
Oe 


an 


| | : 2 tee . 
(g) Use the answer to Question (f) to find an expression for — in terms of —. 
. Cx Cx 


‘(h) Use the answers to Questions (e) and (g) to write down a partial differential equation 
3 A2 AZ 
s C 


wens tee 6T 
involving — and —, |i. eliminate —}. 
Cx Ct Ox 


(i) What can you say about s(0, t)? 

(j) What forces are acting on the mass M? 

(k) What is the total downward force acting on the mass M? 

(1) What is the downward acceleration of the mass M (in terms of the function s)? 


A2 
(m) Use Newton’s second law to find a relationship between T(/,) and =a llo. t). 
a 


(n) What could you say about T(/)) if s were a function of x only? 


~ 


cS. 
(0) What could you say about — when x = /9, if s were a function of x only? 


~ 


CX 


(p) Put s = X(x) in the partial differential equation obtained in Question (h), and hence 
integrate twice to find a solution which is independent of time. 


(q) Do you recognize the solution of Question (p)? 
(Look again at Example 1 of Section 2, page 10.) 


(r) Make a substitution 


a3 (M + Mo)g | Mogx* 
y = s(x,t) — re + i} + “a 


in the partial differential equation obtained in Question (h). What equation do you 
obtain? 


(s) What can you say about y(0,t)? 


(t) Assume now that M = 0. (The case M ¥ 0 is too difficult for us at this stage.) What 
can you say about T(/,)), using the answer to Question (m)? 


A 


C 
(u) What can you say about ie using the answers to Questions (f) and (t)? 
Cx 


A 


(v) What can you say about a using the formula for y in Question (r), and the 
Cx 


} 
4 


answer to Question (u)? 


(w) Use the method of Subsection 3.1 (considering three cases exactly as before) to find a 
set of solutions of the wave equation obtained in Question (r) with boundary 
conditions 


10,1) = 0 
for all f, 


aire 
Cx 


these being the boundary conditions derived as answers to Questions (s) and (v). 
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Appendix: Solutions to the exercises 


Solutions to the exercises in Section 1 


1(i) F(x,y) =~, so 
y 


ie 6 oe 
ox ae eee 
(ii) F(x, y) = Ax? + 2Bxy + Cy’, so 


OF 
— = 2Ax + 2By, 
Ox 


OF 

—— = 2Bx + 2Cy. 

oy 

(iii) F(x, y) = sin(2x + y), so 
OF OF 
— =2cos(2x + y), —— = cos(2x + y). 
Ox oy 


(iv) F(x, y) = o(2x + y), so 


OF 
Fees 2'(2x + y) (where #’ denotes the 
. derived function of @), 
OF 
Wie p'(2x + y). 
y 
(v) F(x,y) = o(ax? + by”), so 
F 
OO = tacdiiax? by); > — Aya ad + by?). 
Ox oy 
2.1) f'(x) 
0 
(ii) ae (x, t) 


a2 


(iii) (x, t) 


OF OF 
3.) F(x, y) = xy, so — = y, — = x, and thus 
ox oy 
ig ee ge pee 
pe ee 
ae 
Oxdy oye 
(ii) F(x,y) =x? + 2x’y + y’, so 
OF OF 
— = 3x? + 4xy, — = 2x? + 3y?, and thus 
Ox oy 
oS a 06°F 
—s = 6x + 4y, = 4x, 
Ox? = Cyex 
07°F 06°F 
Oxéy oy 
(iii) F(x, y) = e**”, so 
OF CF 
— = et’ —_ = e**”, and thus 
Ox oy 


OF oe Or OF 
Ox? @ydx «= Gxdy Gy? 


Indeed, all the partial derivatives of F, of all orders, are equal 
lo“ *”. 


+y 


4(i) U =f(x)g(y), so 


dU dU ee 
—=f'(x)g), =f oe’), —— =S'e)e'0) 
ox oy Oxdy 
oU 0U ; 
Thus, ———— = f(x)g(v)f'(x)g') 
Ox oy 
= 67U 
gay: 
(ii) U =f (ay — bx), so 
U 
i = —bf’(ay — bx) and ee = af (ay — bx). 
ox oy 
oU oU 
Thus, ae 3 + "3 = (), 
(ili) y = f(x —t)+ g(x +f), so 
0 
D = f(x — t+ ex +0) 
Ox 
oy : 3 RC 
a J 1x —1)-+.21(X.+ 0) giving 
o2 
ae =f"(x-—t)+g"(x+t) and 
Ox 
0? 
More nseurd 
Ot : 
ee 
Thus, a —— res 


5. (i) and (iv) are homogeneous linear. (ii) fails to be linear | 

OF OF cee 
because of the term xa oy’ in which two partial derivatives 
are multiplied together, while (iii) fails to be homogeneous 
because of the term e*”, which does not involve the unknown 
or any of its derivatives. The only one of the four which is 
constant-coefficient is (iv). 


Solution to the exercise in Section 2 
1. Substituting s = U + x into Equation (13) gives 


a +x)= 6 + x). 
Thus 

2b Lead Els 

Oe, a ae ee a 
Now, = = 1, so a3 = 0. Also ~ = 0 since ~ means the 
rate of change as t varies and x remains fixed. Thus f= = 0, 


ot? 
and we obtain Equation (14). 


Solutions to the exercises in Section 3 


1.) Substituting c = 2, » = 1 gives 
-T"(t) = 4T(t) and X"(x) = X(x). 
From Unit 6, Section 1, we have 
T = Ae*' + Be~** for some constants A and B;: 
also 
X =Ce*+De™* for some constants C and D. 
(ii) Substituting c = 2, u = —1 gives 
T’=-4T and -X”"=-X. 
Hence 
T = Acos2t + Bsin 2t 
and X =Ccosx + Dsinx 
for some constants A, B, C and D. 
(iii) Substituting c= —2, u = 1 gives 
T" =4T and X”" = X, 


the same equations as in part (i). Thus the solution is the 
same as in part (1). This shows that although the sign of wu is 
important, the sign of c is irrelevant. 


2. We know that Ae” + Be~™ = 0 for all values of t, so 
that, in particular, it is true for t = 0 and t = 1/h. Therefore 


A+B=0 
B e 
Ae+—=0 
e 
and this pair of equations has the unique solution 
A=B=0. 
x Mae | 8 


0(x, t) = X(x)T(t), 


then differentiating partially with respect to x (for fixed t), we 


obtain 
06 
Ox? 

and differentiating partially with respect to t (for fixed x), 


= X"(x)T(t) 


a = X(x)T'(O). 


It follows that 

X"(x)T(t) = c?X(x)T(t), 
so that 

mw 1k) FG 

2) te 


Using the same argument as before we deduce that y is 
constant, and hence the functions X and T satisfy 


A =e ek |. ane. Tut 


4. Using the formula for A,, we have 


2 TT 
= | (x) sin nx dx. 
™ Jo 


Using the hint (with / set equal to z), we obtain 


8d nT 
A. = 
nnn? 2° 


MST204 32 Solutions 


Also B, = 0 for all n, since w(x) = 0, and it follows that 


sin nx cos nct 


jyntl 
rr Ga iP sin (2m — 1)x cos (2m — 1)ect. 


(This is the series used to model the motion of a vibrating 
string in the television programme ‘Waves — Partial 
differential equations’ which accompanies this text.) 


Solution to the exercise in Section 4 


1.4i) Using the boundary condition 
U(L, t) = Rsinat, 


we have 
X(L)T(t) = Rsin at, so that 
R sin wt 
T(t) = 1 
(t) X(L) (1) 


(ii) Substituting the above expression for T(t) into the 
equation T”(t) = c?uT, we have 


Row’ sinwt — c*uRsinat 
X@) 3) 
2 
so that uw = — 
c 


2 
eee ° 63) . . 
(ii) Putting » = ——, in the equation X” = yw X, we have 
c 


2 


aw 
Es paar A — X, 
c2 
so that 
X(x) = Ceos—— + Dsin—. (2) 
‘e 


(iv) If U(O,t) =0 for all t, then X(0) 
Equation (2) we deduce that C = 0. Thus we have 


= 0, and using 


X(x) = Dsin = (3) 
(v) U(x,t) =X(x)T(t) 
= Dsin — Rin wt (using Equations 
Pence (1) and (3)) 
: 3 


<= RA; 
Rsin—sin wt 
Cc 


a 
sin —— 
c 


Solutions to the problems and questions in 
Section 5 


The solutions to the problems are set out below as answers 
to the sequence of questions. The answers to Questions (a) to 
(h) constitute the solution to Problem 1; from (i) to (p) 
constitutes the solution to Problem 2; and from (q) to (w) 
forms the solution to Problem 3. 
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1.(a) The element os was originally of length 6x, so that its 


mass is 
0 


(b)(i) The tension upwards at P; 
(ii) the tension downwards at-Q; 
(iii) the weight of the element 6s acting downwards. 


M 0BOx 


(c) T(x + 6x) — ‘s 


i 52 So ecuee 
(d)(i) 


5 SOs Os 
(ii) a, t) [abbreviated to > 


S(x, t). 


0? 0? 
(iii) 7 (x, t) abbreviated to a I, 
Mogo M,0x 6? 
(e) T(x + 0x) — T(x) + a ee, 
lo ig. 0F 
and in the limit, after dividing by ox, 
OT r Mog 2 Mo és 
a ee ee 
) T- a s i. 


(Use Equation (3) on page 10 with 7) = 0, since the spring is 
originally unstretched; and use a partial derivative because s 
is a function of x and t.) 


OT 6s 

—=k—. 

(g) Ox Ox? 
Os Moe —M, 0’s 
ERS 0 a ied Sell 
(h) Ox? lo lo 6t? 


This is the solution to Problem 1. 


2.(i) The top of the spring is fixed, so that s(0,t) = 0 for 
all t. 


(j) Its weight Mg downwards and the tension in the spring 
upwards. 


(k) Mg — T(Io). 
2 


(i) ; : apz lost in other words — nS 


cas evaluated at (Io, ft). 


2 


rs) 
(m) Mg — T(lo) = M = (lo,!). 


2 


4 would be zero. 


Mg, since in this case —- a2 


(n) T (lo) = 


d 
(o) ~ can be replaced by oa if s is a function of x only, 


and from the answers to Questions (f) and (n) we have 
ds 
Mg =k i— - 1 
P x | 
when x = Io. 


(p) The equation of Question (h) becomes 


k d2X us M o£ 
ee 
so that 
dX Mogx 
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Putting x = /) and using the answer to Question (0) we have 
Mg+k= —Mog + A, 

so that 
A=(M+Mo)gt+k 


and therefore 


pax M 
— = (M + Mo)g + k - OF. 
0 
Integrating again, 
Mogx* 
kX = ((M + Mo)g + k)x — eee B, 


and from the answer to Question (i) we have B = 0, so that 


y= ee i fe } 


This is the solution to Problem 2. 


Mogx? 
eS ie 


3.(q) This is precisely the solution which we obtained in 
Example 1 of Section 2 (see Equation (8) on page 12). In this 
context it is just what we should expect—namely the 
solution which is independent of time. 


(r) The wave equation, 


oe 1 ey 
dx? c? Gt?’ 
k 
where c? = ea 


0 


This is the hint suggested at the end of Problem 3. 
(s) y(0,t) = s(0,t) = 0 for all ¢. 
(t) T(lo) = 0 because M = 0. 


This is the solution to Problem 3(i). 


(u) From the answers to Questions (f) and (t), 
4) 
a = 1 for all t. 
Ox 

This is the solution to Problem 3(ii). 


(v) Using the formula in Question (r) (and putting M = 0), 


oy a Os =e += 
Ox Ox eee 
so that using the answer to aS (u), 
oy Mog Mog 
Plt) =~ F +1)+ k 
asf for all t. 


(w) The discussion of Cases I and II is almost identical to 
that of Subsection 3.1, and there are no non-trivial solutions 
of either form. For Case III we have 


y(x, t) = (Ccoshx + Dsinhx)T(t). 
Since y(0,t) = 0 for all t, we must have C = 0 for a non- 
trivial solution. 


0 
Since = (lo, t) = 0 for all t, we have 


Dhcos hlo = 0, 
so that 
hse 2 + he 


for some non-negative integer n. 
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The required set of solutions is therefore 


lo 


1X 1 act 1 tect 
yn(X, t) = sin(n + 27 A, cos(n + rae + B,sin(n + 5)——}. 
10) 1) 


This is the solution to Problem 3(iii). 


The gerieral solution for y(x, t) is therefore 
ates TX 
Y> sin(n + 4)— 

n=0 lo 

and so, undoing the substitution of Question (r) while 
remembering that now M = 0, the general solution for the 

original function s(x, f) is 

5 Mogx? 


- 1) xe 


t t 
A, cos(n + aes) + B,sin(n + 5 
lo lo 


Mog 
k 


s(x,t) = | 


ae x 
+ >} sin(n + $)— 


T 1 tect 1 tect 
A, cos(n + 4)— + B, sin(n + 3)——}. 
n=0 lo lo lo 


As in Subsection 3.3, a knowledge of the initial conditions 
would allow us to calculate the values of the A,s and B,s, 
thus solving as a Fourier series the problem of a heavy 
spring vibrating vertically under the effect of gravity. 


> r 
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